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ABSTRACT 

We present an investigation of small-scale-structure in the distribution of large 
molecules/dust in the interstellar medium through observations of diffuse interstellar 
bands (DIBs). High signal-to-noise optical spectra were recorded towards the stars 
p Oph A, B, C and DE using the University College London Echelle Spectrograph 
(UCLES) on the Anglo-Australian Telescope. The strengths of some of the DIBs 
are found to differ by about 5 - 9% between the close binary stars p Oph A and B, 
which are separated by a projected distance on the sky of only c. 344 AU. This is the 
first star-system in which such small-scale DIB strength variations have been reported. 
The observed variations are attributed to differences between a combination of carrier 
abundance and the physical conditions present along each sightline. The sightline to- 
wards p Oph C contains relatively dense, molecule-rich material and has the strongest 
/i/15850 and 4726 DIBs. The gas towards DE is more diffuse and is found to exhibit 
weak 'Cy (blue) DIBs and strong yellow/red DIBs. The differences in diffuse band 
strengths between lines of sight are, in some cases, significantly greater in magnitude 
than the corresponding variations among atomic and diatomic species, indicating that 
the DIBs can be sensitive tracers of interstellar cloud conditions. 

Subject headings: ISM: lines and bands — ISM: clouds — ISM: structure — stars: 
individual (p Oph A, B, C and DE) 
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1. Introduction 



Knowledge of the spatial distribution, dynamics, composition and chemistry of neutral and 
ionized gas and dust is of fundamental importance in understanding the interstellar medium (ISM) 
and its relation to star formation. Through observations of absorption lines of atoms and small 
molecules at wavelengths ranging from radio to ultraviolet, the gaseous component of diffuse 
and translucent interstellar clouds has been probed in some detail and indicates that filamentary 
or sheet-like small-scale-structure is common on scales as small as a few AU (Wa tson & Mever 



19961 : lHartquist et al.ll2003c, lHeilesll2007f) . In contrast, the distribution of large molecules and dust 
in the diffuse ISM is far less well characterized, although for denser regions insight has been 
gained through imaging of dust emission. Variations in diffuse i nterstellar band strengths have pre 



viously been reported on distance scales down to a few parsecs (|Points et al.ll2004l : Ivan Loon et al. 
2009k. but oyer di s tance scales ~ 10 — 150 AU, no significant variations have yet been observed 



(ISmith et al.ll2.Q12k iRollinde et al.ll2003|) . The results presented in the present study constitute the 
first detection of DIB variations over distances on the order of a few hundred AU. 



Although the carriers of the DIBs remain to be identified (IHerbigJ 1 1995k ISarrel |2006|) . it is 
widely thought that the absorption features arise from electronic transitions in large gas-phase 
molecules which are associated with, and are possibly formed from, dust grains. The linking of 
diffuse band absorption to the presence of dust grains rests in significant part on the generally good 
correlation between their strengths and E B _ V although deviations are well established. 

In order to probe the large molecule/dust distribution in diffuse regions of Ophiuchus we 
have undertaken observations of diffuse interstellar bands towards four bright stars in the p Oph 
system (A, B, C and DE) - as shown in Figure \T\ - the first results of which were described by 



Cordiner et al.l (|2005L I2006|) . Taking a distance to p Oph A and B of 11 1^ pc (|van Leeuwen 



2007|) . the stars comprising this binary are separated by a sky-p rojected distance of c. 344 AU. The 



stars p Oph C and DE are more distant at 125+[j and 135^}q pc (|van Leeuwenll2007|) . and separated 



from p Oph AB by sky-projected distances ~ 17,000 AU and ~ 19,000 AU, respectively. An 
additional motivation for this work is investigation of the relative band strengths for closely aligned 
lines-of-sight as this could provide a new clue in the search for diffuse band assignments. Of 
particular importance in this regard, the physical and chemical conditions of very closely-separated 
sightlines are expected to be similar to each other, thereby reducing the number of physical and 
chemical variables that may affect diffuse band strengths. In this Letter we describe the detection 
of significant differences in DIB strengths of up to c. 9% between the closely-spaced lines of sight 
towards p Oph A and B, and larger differences up to c. 40% between AB and p Oph C and DE, all 
four of which have the same E B ~v within + 0.01. 
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2. Observations 



The p Oph system consists of five young early-to-mid B-type dwarf stars (|Dommanget & Nys 



19941) . basic properties of which are given in Table [T| Observations were carried out during the 
months of 2004 March and 2005 April using UCLES at the Anglo-Australian Telescope with 
a 1" slit width. Complete wavelength coverage was obtained from 4,600 - 10,000 A, with a 
spectroscopic resolving power of 52, 000 - 58, 000. The median seeing was 0.95" (2004) and 1 .7" 
(2005). This was sufficient for the p Oph AB pair to be resolved on both epochs, with at most 
0.1% cross-contamination of flux in their reduced spectra (assuming Gaussian spatial profiles). 
The separation of p Oph DE is only 0.6", but D dominates the spectrum because it is 1 .2 mag 
brighter than E. 

The search for small-scale-structure through diffuse interstellar band spectra is facilitated 
where there are similarities in stellar spectral type coupled with a low level of contaminating stellar 
photospheric features. The spectral match and stellar spectra are particularly good in the case of 
p Oph A and B, and the situation is generally satisfactory for comparison of AB with p Oph DE. 
The spectrum of p Oph C is contaminated by weak absorption features that are not present for the 
other stars. Short-timescale variability is detected in several of its Sin and He i lines, somewhat 



reminiscent of /3 Cephei-type or slowly-pulsating B stars (see, e.g. iTelting et al.ll2006|) . Contam- 
ination of the observed p Oph C DIBs is negligible however, with the exception of /L15780 and 
5797. 



Reduction and analysis techniques have been described previously in detail by lCordiner et al 



(2006). Particular attention has been paid to non-linearity in the CCD response, scattered light 
subtraction, telluric and blaze corrections. All four program stars were observed in immediate 
succession and without altering the optical setup of the telescope or spectrograph (during each 
night of observations), in order to maximise consistency between the recorded spectra. Exposure 
times were varied so as to obtain similar photon counts for each star; signal-to-noise ratios of the 
reduced spectra are given in Table [TJ 



3. Results and Discussion 

3.1. Spectral observations 

Figure |2] shows the reduced (non-telluric-corrected) spectra (1^ vs. A) in the wavelength region 
covering the /U5780 and 5797 diffuse interstellar bands towards p Oph A, B, C and DE. The lower 
three panels show the ratios: 1 A (B) / 1 A (A), 1 A (C) / 1 A (A) and I A (DE) / 1 A (A). The ^5780 and 5797 
DIBs are stronger, by 5% and 9%, respectively, towards B than A, as indicated by the residuals in 
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the L;(B) / ^(A) plot that precisely match the positions and profiles of the respective DIBs. These 
two bands are also stronger towards DE than A, but in this case /15780 exhibits a much greater 
change in strength of ~ 20%, the change for A5191 being only c. 7%. This is a notable finding 
given that the values of Er-v towards p Oph A and p Oph DE are almost identical with values of 
0.47 and 0.48 ( Pan et al. 20041) . Despite contamination of the p Oph C spectrum, it is apparent that 
both /15780 and A5191 are stronger towards C than A. 

Figure |3] (top panel) contains the corresponding data for /16614 which again show this DIB to 
be slightly stronger towards p Oph B than A, and much stronger towards DE than A, in a similar 
fashion to AA5191 and 5780. A significant decrease in /16614 band strength is seen towards p Oph 
C relative to p Oph A. The /J6614 spectra appear to show some evidence o f profile variation which 
may result from different carrier internal temperatures (ICami et al.l l2004) or isotopic enrichments 
iWebstei]ll99~6l) . 



Measurements of the equivalent widths for AA57S0, 5791, 6614 and seven other diffuse bands 
are given in Table |2] Broader DIBs are omitted due to the difficulty in accurately defining their 
continua in echelle spectra and the increased likelihood of overlapping stellar features. The results 
presented here for the yellow-to-red DIBs (obtained in 2005) match closely the results obtained by 
Cordiner et al.l (|2006|) . who used only the 2004 data (which had slightly lower signal-to-noise, but 
was observed with better seeing), and focussed on DIBs at wavelengths > 5500 A. The fact that 
the spectra obtained on these two epochs match so closely despite the different instrumental setups 
and observing conditions adds confidence to the results presented here. 

The /M5780, 5797, 5850, 6379 and 6614 DIBs have the greatest central depths of those in our 
sample, and are the only ones to show a significant difference in strength between p Oph A and B. 
Small differences were measured for other bands towards these stars, but they are not significant 
given the noise level. Our data are consistent with all DIBs being slightly stronger towards p Oph B 
than A (by c. 5%). However, for the AB pair vs. C and DE, the situation is more complex. Out 
of all the p Oph sightlines, we find that /J5850 is strongest towards C, whereas /16614 is weakest. 
Additionally, DE has clearly the strongest ,16196 DIB, but the weakest /14726 (the latter DIB is 
40% weaker towards DE than C and 32% weaker than AB). This wide variation in behavior of 
different DIBs between the closely-spaced sightlines in our sample shows that the carriers are 
being influenced by quite subtle variations in physical/chemical interstellar conditions. 

An important advance in diffuse interstellar band studies was the identification of a set of DIBs 
that have a propensit y to follow in strength the column density of diatomic carbon molecules - the 
so-called 'C 2 DIBs' (IThorburn et al] 120031) . A spectroscopic property of the C 2 DIBs (possibly 
related to the structure of their carriers), is that they fall predominantly at the blue end of the 
spectrum (< 5500 A), whereas the non-C2 DIBs that we have measured fall at the red end (> 5500 
A). Among the strongest of the C 2 DIBs is AA96A, data for which are shown in Figure |3] (bottom 
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panel). This DIB shows no detectable difference between p Oph A and B but is found to be weaker 
towards p Oph DE than p Oph A. Interestingly, we find that a total of nine out of the ten observed 
C 2 DIBs (/U4726, 4 734, 4964, 49 84, 5176, 5418, 5512, 5541, 5546) are measured to be weaker 
towards DE than AB ( Smith Il2006|) . The remaining /15170 C 2 DIB is overlapped by abroad spectral 
feature, probably stellar in origin, which precludes a proper measurement of its variation. Ten out 
of eleven of the non-C2 DIBs we observed in the yellow-to-red part of the spectrum (for which 
measured variations are larger than the errors), show the opp osite trend, an d are stronger towards 
DE than AB (data for X15705, 6660 and 7224 are given bv ISmith lEoOol) . the single exception 
being /J5850. These results strongly support the association of the C 2 DIBs as a special group. 



3.2. The ISM towards the p Oph star system 



Snow et al.l (120081) studied the atomic (K i, Na i and Ca i) column densities towards p Oph A, C 



and D and determined the three-dimensional locations of these stars with respect to the p Oph cloud 
complex. They conclude that while p Oph AB and D are probably in front of the dense molecular 
cloud and probe a more diffuse portion, p Oph C more likely lies embedded in a CN-rich molecular 
part. As discussed by Snow et al. (2008), the ratio of total-to- selective extinction (R v ) is unusually 
high towards p Oph, and the far-UV rise is quite flat. These observations are cons istent with the 



prese nce of relatively large dust grains and an absence of the very smallest grains (|Seab & Snow 
1995J). Notably, the reported E B _ V values for the p Oph stars are almost identical (Table [T), which 



given the magnitude of the observed DIB variations, is consistent with prev ious finding s that the 
DIBs do not arise directly from the large grains that cause optical extinction (lSarrel l2006). 



3.3. Comparison of DIBs with known atoms and molecules 



A body of data for atoms and small molecules towards the p Oph stars is given in Table 
|2l where the column densities quoted are in general the sum over a number of blended velocity 
components. Unfortunately, the widths of the DIBs preclude their association with any particular 
component. We measured equivalent widths for 23 lines of C 2 in each of our spectra (in the range 
7700 - 8900 A, which covers the (2-0) and (3-0) C 2 vibrational bands) . Total C 2 column densities 
were calculated using the oscillator strengths of iBakker et al.l (|1997|) : the best-fitting rotational 
temperatures (r rot (C 2 )) were obtained using B. McCall's C 2 calculator (at http://dib.uiuc.edu/c2/), 
and are given in Table |2J All four lines of sight show absorption due to CN at around 2 kms" 1 , 
with a common Doppler b ~ 0.& km s" 1 . 



Apart from CH + which has a larger total column density towards p Oph A than B, the errors 
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on the total column densities are too large to permit any clear distinction to be drawn between the 
properties of the gas towards p Oph A and B. However, if one considers only the main 2 kms" 1 
component, CH, CH + , Ki and Can are all 10-20% stronger towards B than A. Our high signal-to- 
noise Ki equivalent width measurements also show more neutral potassium towards B than A. We 
infer that the column of DIB carriers increases towards p Oph B in concert with the species noted 
above. It is difficult, however, to know whether this is due to variations in the gas density, lower 
atomic (and DIB) depletions, changes in ionization equilibrium, or due to the geometric structure 
of the cloud. 



Our observations towards p Oph C, as well as those of lPan et al.1 (I2004|) and lSnow et al.1 ([2008) 
show a relative under- abundance of K i and Cai com pared with t he o ther three sightlines. Consis- 
tent with its (~ 3 times) greater CN column density, ISnow et al.l (|2008|) concluded that p Oph C is 
more deeply embedded in the dense, CN-rich material of the background (star-forming) molecular 
cloud. This relatively high density material should result in more depletion, which is the most likely 
explanation for the low Ki and Cai column densities. The lower /16614 strength towards p Oph C 
may also be related to this depletion. The fact that the variations in the 'red' DIB strengths are 
only on the order of a few percent for C compared with AB and DE shows that their carriers are, in 
general, not closely associated with CN. This result is in accordance with the study of A/15797 and 
5780 by lWeselak et all (|2008l) . By contrast, the 'blue' ,14726 DIB is very st rong towards C whic h 
confirms the association of its carrier with molecular material traced by C2 (IThorburn et al.ll2003|) . 



The greater strengths me asured for AA5797 and 5850 towards p Oph C are consistent with the 
assignment of these DIBs to iKrelowski & Walked (I1987|) 's group III, which is a family of diffuse 
bands that tend to have relatively narrow profiles. These 'KWIII' DIBs (in contrast to those in the 
KWII group) are generally observed to be strongest in better shielded, de nser diffuse clou ds where 
the UV radiation field is weaker and small molecules are more abundant dCami et al. 1997 ). 



The sightline towards p Oph DE is relatively m ore diffuse tha n C, as indicated by the lower CN 
column density. It has a lower depletion factor (F*; Jenkinsll2009 ) than A, indicati ve of a lower ga s 
density. The 2 km s~ ! component towards DE also has the lowest modeled density (|Pan et al. 
consistent with it being the more diffuse sightline. The greater C 2 rotational temperature measured 
towards DE may reflect the increased (photo-electric) heating expected to result from a stronger 
UV radiation field pervading this more diffuse gas. As shown in Figure|21 the amount by which the 
A51S0 DIB strength increases between A and DE (~ 20%), is greater than t he amount by which 
A5191 increases (~ 10 %). This pattern is consistent with the hypothesis (|Krelowski & Walker 



19871 : ICami et alJll997|) that the carrier of the ,15780 DIB favors more diffuse environments than 
,15797. The strengths of other DIBs in the same (KWII) group as A51S0 0M6196, 6203 and 6614) 
show a similarly large enhancement (~ 20%) towards DE compared to A/B. With the striking 
exception of CH + (-26%), CH, Cai, Ki and Can show increases in column density in the range 
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15-30% (for the main 2 kms" 1 component) from A to DE, which is analogous to their increase 
from A to B discussed earlier, and could plausibly be related to a lower degree of depletion in the 
gas towards DE. 

The relative weakness of the C2 DIBs towards DE is consistent with this being the more 
diffuse (and therefore, more molecule-poor) sightline of our sample. Our C 2 spectra are, however, 
of insufficient signal-to-noise to discern differences in the C 2 column densities between the p Oph 
sightlines; dedicated future C 2 observations will be required to confirm this. 



4. Conclusions 

Using high signal-to-noise echelle spectra of the stars p Oph A, B, C and DE, the strengths of 
some of the observed diffuse interstellar bands are found to differ by about 5% on sky-projected 
distance scales less than c. 344 AU. The observed DIB strength variations between p Oph A and 
B are attributed to an increase in the abundances of their carriers/chemical precursors from A to B. 
Towards p Oph C and DE, the DIB behavior is more complex, probably as a result of the greater 
differences in physical and chemical conditions in the ISM over the relatively larger distances 
between these sightlines and the different locations of these stars with respect to the background 
molecular cloud. The yellow/red DIBs are generally found to be strongest in the relatively diffuse, 
less-depleted gas towards p Oph DE. This is in striking contrast to the C 2 (blue) DIBs, which are 
consistently weaker towards DE than the other stars. The /15850 DIB is found to be strongest 
in the more d ense, depleted and molecul e-rich diffuse gas towards p Oph C, consistent with its 



assignment to lKrelowski & Walked (|1987|) 's group III. 



Observed DIB strength variations are, in some cases, significantly greater in magnitude than 
the corresponding variations among atomic and (most) diatomic species; for example the equiva- 
lent widths of /M4726, 5780 and 6614 towards DE compared with A/B exceed the magnitude of 
variation in all observed species apart from CH + . Such variations imply that diffuse band carrier 
abundances are very sensitive to variations in cloud conditions. Given the small distances between 
our observed sightlines, differences between the elemental compositions of their gases are expected 
to be negligible. Variations in diffuse band strengths and atomic and molecular column densities 
could therefore arise as a consequence of density inhomogeneity or the strength of the incident 
radiation field. Recent observations of CH + and SH + in the diffuse ISM have highlighted the im - 
portance of turbulent dissipation arising from shocks or velocity shears (bodardetall2009[l2012l) : 
this may also play a role in determining diffuse band carrier abundances on small distance scales, 
particularly if the carrier formation is driven by warm chemistry or through release of material 
from interstellar dust grains. 
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A more complete explanation for the origin of the observed small-scale variations in DIB 
strengths in the p Oph system will require a detailed knowledge of the properties of each sight- 
line, including the density, temperature, (gas and solid-phase) chemical abundances and radiation 
field. This could be achieved through high signal-to-noise, high-resolution optical/UV/IR studies 
of atomic and molecular absorption lines including Hi, H 2 , C, C2, and other neutral and ionic 
species sensitive to physical conditions in the ISM. 

We thank PATT for the allocation of AAT time and T & S, EPSRC for studentships, and STFC 
for a visitor grant. MAC thanks the NASA Astrobiology Institute through The Goddard Center for 
Astrobiology. SJF thanks M. M. Dworetsky and I. D. Howarth for discussions on photospheric 
line-profile variations in p Oph C. 
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Fig. 1 . — Photographic J (blue) image of the p Oph system obtained from the Digitized Sky Survey 
(DSS). Angular distances between components are as follows: #(AB) = 3.1", #(AC) = 150", 
0(AD) = 169", 0(DE) = 0.6". 
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Table 1 . Program star properties 



Star 


HD 


Sp. Type 


V 


Eb-v 


S/N (6000 A) 


p Oph A 


147933 


B2/B3 V 


5.0 


0.47 


1200 


pOphB 


147934 


B2 V 


5.9 


0.47 


1100 


pOphC 


147932 


B5 V 


7.3 


0.47 


800 


pOphD 


147888 


B3/B4 V 


6.8 


0.48 


900 



Note. — Spectr al types are from the Michig an Catalogue 
of Spectral Types (Houk & Smi th-Moorel Il988|) . Reddening 
values {E B -v) are from lPan et al.l (|2004[) . 
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Table 2. p Oph diffuse interstellar band and atomic/molecular line measurements 





A 


B 


C 


DE 


Error 


Equivalent widths (mA) 


/14726 


99.3 


97.8 


112.9 


67.4 


4.3 


/14964 


23.5 


22.0 


21.5 


19.5 


0.8 


,15780 


188 


197 


[199] 


227 


2.8 


A5191 


49.7 


54.2 


[56.0] 


55.2 


1.8 


,15850 


29.8 


31.9 


35.2 


25.2 


1.8 


,16196 


15.9 


15.6 


15.0 


18.8 


2.0 


,16203 


42.0 


42.5 


43.0 


47.2 


6.0 


/16376 


14.8 


15.8 


15.5 


21.0 


1.3 


,16379 


27.1 


28.3 


28.5 


34.7 


1.3 


,16614 


63.9 


67.7 


62.1 


84.9 


1.3 


K 1/17699 


88.1 


89.5 


77.7 


90.5 


0.3 


Column densities and derived values 


Af(Hi)/10 21 


4.3 






4.9 


1 


N(H 2 ) / 10 2() 


3.7 






3.0 


1 


Af(Ki)/10 12 


1.03 


1.09 


0.80 


1.10 


0.05 


N(Cai)/ 10 10 


1.74 


1.69 


1.62 


2.02 


0.09 


^V(Can)/ 10 12 


1.75 


1.96 


1.92 


1.94 


0.1 


N(CU) / 10 13 


2.37 


2.33 


1.92 


2.18 


0.1 


N(CU + ) / 10 13 


1.59 


1.43 


0.70 


0.75 


0.06 


N(CN) / 10 13 


0.21 


0.20 


0.60 


0.21 


0.01 


N(C 2 ) / 10 13 


4.1 


2.9 


3.3 


4.0 


1 


r rot (C 2 )/K 


51 


49 


46 


56 


1 


n 1 cm -3 


625 


450 


1100 


425 




Depletion (F*) 


1.09 






0.88 


0.07 



Note. — Upper part of table gives equivalent widths of inter- 
stellar features measured towards p Oph stars; values in brack- 
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ets indicate probable stellar photospheric line contamination. 
Lower part of table gives atomic and molec ular column densi- 
ties (in cm -2 ) taken from Pan et al.1 (|2004|) except: C2 (from 



this w ork), H 1, H 2 and depletion factor (F*) from Ijenkins 
( 2009 ). Gas number densities n are from the chemical model 



-1 



of iPan et al.l (|2005|) and pertain to the CN-containing 2 km s 
component only. Final column gives an average (+) error esti- 
mate on the respective values. 
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5770 5780 5790 5800 5810 

Wavelength (A) 



Fig. 2. — Spectra of the region covering the /L15780 and 5797 diffuse bands. Top four panels 
are observations towards p Oph A, B, C and DE, with ordinate offsets of 0.1; these are corrected 
for the blaze function but not for telluric effects and are unsmoothed. The three lower panels are 
non-blaze-corrected non-telluric-corrected ratios I /1 (B)/I /1 (A), I /l (C)/I /l (A) and l A (DE)/l A (A). The 
spectrum of p Oph C is significantly contaminated by weak stellar lines in this region. For a ratio 
X/A, a value < 1 means the band is stronger along sightline X than A. 
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6600 6610 6620 

Wavelength (A) 




4955 4960 4965 

Wavelength (A) 



Fig. 3. — Spectra of the region covering the /16614 diffuse band (top) and /14964 'C2' diffuse band 
(bottom). Data are as described in Figure 0but with ordinate offsets of 0.05. 



